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I.  INTRODUCTION 


The  acoustic  environment  of  an  engine  test  cell  of  the  type  at  the 
Engine  Test  Facility  of  AEDC  is  extremely  complex.  x The  sound  field  present 
in  such  a cell  consists  of  a direct  component,  due  to  the  sound  pressure  waves 
emanating  directly  from  the  distributed  source,  and  a reverberant  field  assoc- 
iated with  waves  reflected  from  the  surfaces  of  the  enclosure.  The  relative 
contribution  of  these  two  components  must  be  determined  in  order  to  either 
model  the  acoustic  field  in  the  test  cell  or  infer  the  strength  and  nature  of 
sources  from  in-cell  measurements  to  predict  the  free  field  sound  that  such 
sources  would  produce. 

One  factor  that  affects  the  relative  contribution  of  the  direct  and 
indirect  fields  is  the  strength  of  the  near  field  component  of  the  direct 
radiation  field.  It  is  well  known  that  the  sound  pressure  near  an  acoustic 
source  can  vary  much  more  strongly  than  the  1/r  relation  that  applies  in  the 
far  field.  In  this  region  the  sound  pressure  is  not  in  phase  with  the  acoustic 
particle  velocity  and  cannot  be  directly  related  to  acoustic  energy  flux.  Thus, 
one  of  the  problems  that  must  be  addressed  before  substantial  progress  can  be 
made  toward  an  understanding  of  the  acoustic  environment  in  an  engine  test  cell 
is  that  of  interpreting  near  field  acoustic  data  and  their  relation  to  the 
acoustic  sources  and  far  field  sound. 

The  primary  objective  of  the  present  study  is  to  determine  the  relation- 
ship between  acoustic  measurements  in  the  near  field  and  the  far  field  sound 
produced  by  a free  circular  jet.  A longer  range  objective  of  the  research 
effort  is  to  attempt  to  establish  techniques  for  predicting  the  far  sound  field 
using  in-cell  acoustic  measurements. 

The  sound  source  used  throughout  this  study  was  a free  turbulent  jet 
exhausting  from  a convergent  circular  nozzle.  In  the  engine  test  cells  at  AEDC, 
the  exhaust  jet  from  the  engine  is  captured  by  a scavenging  duct  or  diffuser. 

The  inlet  of  this  duct  is  normally  located  a very  short  distance  downstream  of 
the  engine  exhaust  nozzle.  In  this  investigation  no  attempt  was  made  to  take  • 
into  account  the  effects  of  the  scavenger  duct  of  the  surrounding  test  section 
walls. 


The  general  approach  taken  in  this  study  was  to  model  the  jet  source 
region,  for  a given  frequency,  with  a combination  of  uncorrelated  elementary 
sources  derived  from  the  multipole  expansion.  The  strengths  of  the  various 
source  components  at  a particular  frequency  were  determined  from  measured  1/3 
octave  far  field  directivity  patterns.  Initially,  it  was  assumed  that  the  jet 
noise  source  distribution  could  be  replaced  by  an  equivalent  stationary  source 
distribution  which  could  be  approximated  by  a group  of  elementary  sources 
centered  in  the  nozzle  exit  plane.  During  the  course  of  the  investigation,  it 
was  found  that  in  order  to  closely  approximate  the  near  acoustic  field,  it  was 
necessary  to  introduce  the  motion  of  the  sources  and  their  distribution  along 
the  jet  axis  with  frequency. 
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II.  STATIONARY  SOURCE  MODEL 


In  this  section  equations  are  developed  for  synthesizing  the  source 
region  for  an  axisymmetric  turbulent  jet,  assuming  that  the  true  noise  sources 
can  be  replaced  by  an  equivalent  stationary  source  distribution. 


2.1  MULTIPOLE  EXPANSION 


The  sound  pressure  field  produced  by  a stationary  distributed  source 
f (x^XjjX-jjt)  is  governed  by  the  non-homo geneous  wave  equation 

M^  + M)-  + My ^ -Mr  - - f(x  ,k  x t)  (2.1) 

3x^  3x  a 3t  ±11 

12  3 0 

In  the  absence  of  any  reflecting  walls  the  free-field  solution  to  this  equation 
is  given  by 


p(x,t)  = 


I 


f(y. 


t - R/a  ) 
o 


47X 


dy 


(2.2) 


where  x is  the  position  vector  to  the  field  point,  y is  the  position  vector 
to  a source  point  and  R = - ^| . The  geometry  under  consideration  is 

illustrated  in  Fig.  1 showing  the  various  position  vectors.  It  is  assumed 
that  the  distributed  source  is  located  in  the  neighborhood  of  the  origin  of 
the  coordinate  system. 

The  integral  solution  (2.2)  is  extremely  difficult  to  evaluate  even  if 
the  source  function  is  known  a priori.  However,  by  expanding  the  integrand 
in  a Taylor  series  about  it  = i,  a series  expression  for  the  sound  pressure 
field  can  be  obtained  which  is  frequently  more  tractable.  (Ref.  1 ). 


p(x,t)  = 


z 

jk2  = 0 


,j+k+2 


. j+k+  2 


3J (-1)-  x 


4.x  <2-3> 

J o 


This  series  is  called  a multipole  expansion  and  functions  m^^(t)  are  called 
multipole  moments  of  the  source  distribution: 

W0  = / jj  ki  Ii  f<y’°  dy 


/ 


(2.4) 


8 


AEOC-TR-79-36 


The  multipole  expansion,  in  effect,  replaces  the  distributed  source  by  an 
infinite  combination  of  elementary  point  sources  located  at  the  origin  of 
the  coordinate  system.  The  leading  term  of  Eq.  2.3  where  j = k = l = 0 
represents  a monopole  source  with  sound  pressure  field 


1 

4irx 


mooo 


(2.5) 


The  next  three  terms  of  the  series  with  j+k+  & = 1 correspond  to  dipole  sources 
with  sound  pressure  fields 


3 


The  six  terms  corresponding  to  j+k+i 
ten  terms  with  j+k+f,  = 3 are  known  a, 
j+k+£  = n correspond  to  multipoles  o 


= 2 are  known  as  quadrupole  terms  and  the 
i octupole  terms.  In  general,  terms  with 
: order  2n . 


If  it  is  assumed,  for  the  moment,  that  all  parts  of  source  region  are 
oscillating  in  phase,  at  a single  frequency,  the  source  function  can  be  written 
as 


f (y,t)  = F (y)  a** 


(2.7) 


and  the  multipole  moments  become 


where 


mjk£(t)  ” Mjk£  e 


iwt 


V j A y2  A d? 
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Introducing  spherical  polar  coordinates  as  shown  in  Fig.  2 with 


— = cosGsin^1 


x2  X3 

— = sin6sin$  — = cos  <J 

x x 


the  sound  pressure  fields  for  the  nonopole  and  dipole  terms  can  be  written  as 


n M000  l(Wt  " a ) 

pmu’t)=-^re 


(2.8) 


P,  (x,t) 
1 


i k sin(j)cos0M  . i(ojt — — ) 

XUU  1 v 3n 

4 vx  kx  e 


PH  (*,t) 
d2 


i k sin<f'9inGM-,  . . i(o)t — — ) 

5i0  (1_  i j a 

4ttx  kx 


(2.9) 


i k cos  <Js  M , i(£ut-  — ) 

UU1  / a lx  Ba 

4^ (1“  k^e  0 


where  k = — = -y-  is  called  the  wave  number, 
a A 
o 

It  should  be  noted  that  the  factor  (1  - i/kx)  in  Eq.  2.9  contains  the 
near  field  contribution  for  the  dipole  sound  pressure.  As  kx  = 2’rx/A  becomes 
much  larger  than  unity,  so  that  the  observer  is  many  wave  lengths  from  the 
source,  this  factor  approaches  one,  resulting  in  the  far  field  expression  for 
the  dipole  sound  pressure.  Equations  for  the  quadrupole  and  octupole  contri- 
butions to  the  sound  field*,  similar  to  Eq.  2.9,  have  been  derived  but  these 
are  much  more  lengthy  and  are  therefore  not  presented  in  the  body  of  this 
report.  A full  listing  of  the  equations  used  in  computing  the  near  and  far 
sound  fields  is  presented  in  Appendix  A. 


2.2  COMBINATION  OF  UN CORRELATED  SOURCES 

Considering  the  source  region  as  the  axi-syimnetric  turbulent  jet  flow 
from  a circular  nozzle,  the  sound  field  produced  by  this  source  distribution 
should  possess  axial  symmetry.  As  will  be  shown,  however,  care  must  be  taken 
as  to  how  this  symmetry  condition  is  introduced.  If,  for  example,  the  source 
function  given  in  Eq.  2.7  is  assumed  to  be  symmetric  about  the  x^  axis,  the  axis 
of  the  jet,  it  takes  the  form 


f{y,t)  = F(r,y3)e^t 


(2.10) 
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where 


yl+  y2 


If  such  a source  distribution  is  used,  many  of  the  raultipole  moments  vanish 
identically.  In  particular,  the  dipole  moments  M^gg  and  MqiO  are  zero  and  the 
dipole  contribution  to  the  sound  pressure  field  becomes 


P 


d 


I k COSfjMgg^ 
4irx 


This  pressure  field  would  be  added  to  that  for  the  monopole  and  surviving  higher 
order  poles  to  produce  a sound  field  that  is  axially  symmetric. 

Further  reflection,  however,  reveals  that  the  assumptions  made  in  writing 
Eq.  2.9  are  much  too  restrictive.  At  a particular  frequency,  both  the  source 
region  and  the  sound  field  produced  would  be  axially  symmetric  on  an  instantaneous 
basis  with  respect  to  amplitude  and  phase.  The  acoustic  sources  In  a jet  flow 
are  directly  related  to  flow  turbulence  and  are  therefore  random  in  nature. 

Thus  while  the  acoustic  source  distribution  at  a particular  frequency  may  be 
regarded  as  axi— symmetric  with  respect  to  average  amplitude,  it  is  certainly 
not  symmetric  on  an  instantaneous  basis  as  was  assumed  in  Eq.  2.9, 

A second  method  that  appears  more  realistic  is  to  model  the  source  with 
an  assemblage  of  uncorrelated  elementary  sources  (monopole,  dipole,  etc.)  These 
sources  can  be  combined  in  such  a way  as  to  produce  a sound  field  that  is  axi- 
symmetric  with  respect  to  mean  square  pressure.  This  technique  has  been  used 
by  Mani  et.  al.  (Ref.  2 ) to  model  the  source  distribution  for  a small  circular 
jet  in  an  open  jet  wind  tunnel. 

To  illustrate  this  approach,  consider  the  contribution  to  the  sound  field 
due  to  dipole  radiation.  If  it  Is  assumed  that  the  three  dipoles  given  by 
Eqs.  2.8  are  uncorrelated,  the  mean  squares  of  the  acoustic  pressures  add  to 
give  the  mean  square  pressure  due  to  dipole  radiation: 


(M^qq  sin2  * cos2  9 + M2^sin2  $ sin20 

+ M2Q1  cos2*)  (2.11) 


2 2,2,2 
p = p,  + p.  + p, 
d *'q1  ^d2  ^ 


(4ux)' 


1 + 


(kx)2 J 
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The  sound  field  given  by  this  equation  is  not  symmetric  about  the  jet  axis  for 
arbitrary  values  of  the  dipole  moments.  However,  if  the  strengths  of  the 
and  *2  dipoles  are  the  same  so  that  = ^010*  t*ie  ^ePendence  on  azimuthal 

angle  0 disappears  yielding 

~1  kZ  fl  1 1 

pd  = 1 1 

(4ti  n)  L (kx) 

This  is  the  most  general  form  for  a combination  of  uncorrelated  dipoles  that 
produces  an  axially  symmetric  sound  field.  In  the  acoustic  far  field,  where 
kx  >>  1,  Eq.  2.12  reduces  to 


(Mioos±n2  ♦ + “ooi  cos2 


(2.12) 


(4ttx) 


(M100  3i”2  ♦ 


ajd!  COS2#) 


(2.13) 


Equations  similar  to  Eq.  2.12  for  the  mean  square  sound  pressure  for 
uncorrelated  quadrupoles  and; octupoles  are  given  in  Appendix  A.  The  far  field 
approximations  for  the  mean  square  sound  pressure  for  monopole,  dipole,  quadru- 
pole,  and  octupole  sources  are  shown  below. 


-s-  M 


000 


(4lTx) 


T2 


-tt  , 2 


(4xx) 


2 tM2100sin2  $ + M2q01  cos2*] 


(2.14) 


-r  ,4 


(4rx) 


2 4 2 2 2 2 4 

j [M  20Qsin  <t>  + M llQ  sin  9 cos  q + X 0Q2cos  4»1 


k6  [M2030sin6  (J>  + MZ021  sin4  * cos2  <j>  + M^^sin2  j>  cos4  9 + 


(4irx) 


r,2  6. . 

M 003cos  *] 


The  total  mean  square  pressure  is  obtained  by  combining  the  contributions 
of  the  above  sources,  assuming  that  they  are  uncorrelated.  Thus, 


2 


P 


(2.13) 
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2.3  UNIQUENESS  CONDITION 


In  order  to  complete  the  synthesis  of  the  acoustic  sources,  the 
multipole  moments  M2q00'  m2100’  m2001»  eCc' » Tnus.t  be  determined  for  a given 
frequency.  This  was  done  by  using  a least  squares  technique  to  determine  the 
constants  that  best  fit  the  far  field  sound  pressure  data  at  a particular 
frequency.  The  multipole  moments,  however,  are  not  uniquely  determined  by 
this  procedure.  To  illustrate  this,  suppose  the  far  field  sound  at  some 
frequency  is  accurately  represented  by  the  combination  of  a monopole  and  two 
dipoles.  Then 


(2.16) 


2 2 

By  using  the  identity  sin  $ + cos  4 


1,  Eq.  2.16  can  be  rewritten  as 


2 it  *2  2 *2  2 . 

p = ——2  IH  1Q0  sin  $+  M Q01  cos  *] 

(4irx) 

where 

M2 

M*2  = m2  + %00 

100  100  ,2 


(2.17) 


M*2  = m2  + %00 

001  001  ,2 


2 2 

Similarly,  if  M > M ^qq>  tbe  mean  square  pressure  can  be  written  as 


1 2 

— M 4 2 

2 000  „ K „'2  2± 

P "W  K °“"5* 


<2.18) 


where 


* 2 2 2 
M = M - M 

001  001  100 


m’2ooo  ■ ^ooo  + k2H2ioo 


While  these  three  expressions  give  the  same  sound  pressure  in  the  far 
field,  the  near  field  pressures  are  not  the  same,  because  the  near  field  terms 
are  different  for  a monopole  and  a dipole. 
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In  order  to  uniquely  determine  the  multipole moments,  an  additional 
condition  suggested  by  Man!  (Ref.  2 ) is  imposed  based  on  the  intuitive 
concept  that  the  lower  order  sources  should  have  their  maximum  possible 
strength.  Using  this  condition  in  the  example  above  specifies  Eq.  2.18, 
since  M'qqq  is  the  largest  passible  monopole  moment. 

As  a further  example  of  the  application  of  this  uniqueness  condition, 
suppose  that  the  far  field  sound  is  accurately  expressed  in  terms  of  a 
combination  of  uncorrelated  quadrupoles.  Then 


2 4 2 2 A 

p = Q^sin  <(>  + Q^sin  4 cos  $ + cos  $ 


(2.19) 


where 


k M200  n = R Ml10 

Q1  ,2  ■ Q2 


kV 


(4irx) 


(4ttx) 


2 ’ Q3  ” 


002 


(4ttx)‘ 


4 2 2 4 

The  identity  sin  4 + 2 sin  4cos  4>  + cos  $ = 1 can  be  used  to  extract  the 
largest  embedded  monopole.  Let 


M 


000 


P1  “ 111111  ^1*  “ 471  x 


and  redefine 


Ql  = Ql  - Pr  Q2  = Q2  ' 2?!*  Q3  = Q3  - pi 


Then  Eq.  2.19  becomes 


2 ^ 2 2 ^ 
p = + Q'^sin  <j>  + Q'^sin  4>  cos  + Q'^cos  :{> 


(2.20) 


2 2 

Next  use  the  identity  sin  $ + cos  <£  = 1 to  remove  the  embedded  dipoles. 


Let 


kV 


D1  = min  (Q'r  q'2)  = 


100 


(4vrx) 


and  redefine 


Tt  t 


"l  ' Q1  - Dl’  **2  ' % ‘ V ^ " Q3 
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Similarly  let 
D. 


k2M^ 


2 nin  (Q^  qp 


001 


(4ttx)‘ 


and  redefine 


This  gives  the  desired  representation  for  the  sound  pressure  field  with  the 
lower  order  sources  having  their  maximum  strength 


p2  = + D^sin2$  + D^cos2!);  + 

ii » 4 hi.  2 2 4 

Q.  sin  $ + Q,  sin  <|>cos  <p  + cos 


(2.21) 


This  procedure  is  easily  extended  to  any  level  of  the  multipole  expansion. 
In  the  present  study  the  expansion  was  truncated  at  the  octupole  level,  so  that 
a least  squares  fit  was  used  to  determine  constants 


kV 


030 


(4irx) 


k6M2 


021 


(4vx)i 


v6*,2 
k M 


A3  = 


012 


(4ttx)j 


k6M2 


\ = 


003 


(4ttx)' 


in  the  equation. 


p2  = A^sin^  + A2sin^  $ cos2(f>  + A^sin^  cosA<|> 

+ A^cos^ 


(2.22) 


15 


AEDC-TR-79-36 


The  identities 


6 4 2 2 4 6 

sin  $ + 3 sin  cos  $ + 3 sin  16  cos  ^ + cos  <f>  = 1 


sin^$  + 2 sin^iicos^ij;  + cos^$  = 1 


2 2 
sin  ({i  + cos  <p  = 1 


were  used  to  extract  the  embedded  monopole,  dipoles,  and  quadrupoles  using 
the  procedure  discussed  above. 


2.4  CONSTRAINED  LEAST  SQUARES  PIT 

It  should  be  noted  that  the  constants,  A^,  A^,  A_,  in  Eq . 2.22  must 
all  be  non-negative  according  to  the  derivation.  Consequently,  a simple 
least  squares  procedure  cannot  be  used  to  determine  these  quantities . In 
(Ref.  3 ) , a constrained  least  squares  algorithm  based  on  the  Kuhn  Tucker 
Theorem  is  presented.  A Fortran  subroutine  based  on  this  algorithm  was 
adapted  to  the  present  problem  to  determine  values  of  A^,  k^,  A^»  that 
best  fit  the  measured  sound  pressure  data  subject  to  the  constraint  A^  >_  0. 
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III.  EXPERIMENTAL  PROGRAM 


3.1  EXPERIMENTAL  FACILITIES 

All  acoustic  data  obtained  in  this  study  were  gathered  in  the  UTSI 
free  field  facility  shown  in  Fig.  3.  This  facility  has  been  described  in 
detail  in  Ref.  4 and  only  the  essential  features  will  be  discussed  here. 

The  test  facility  was  basically  designed  to  be  used  to  investigate  the  sound 
field  generated  by  high  velocity  turbulent  air  jets.  The  air  supply  system 
for  the  test  rig  is  capable  of  delivering  air  mass  flows  in  excess  of 
40  lb  /sec  at  stagnation  temperatures  up  to  lOOO'F.  Being  a blow-down 
facility,  the  runtime  is  inversely  proportional  to  the  air  mass  flow  rate, 
with  a flow  of  40  lb^/sec  corresponding  to  about  a 3 minute  run. 

Acoustic  data  on  air  jets  are  normally  taken  with  the  stilling  chamber 
and  exhaust  jet  axis  oriented  at  an  angle  of  about  35  1/2°  to  the  ground 
plane,  as  shown  in  Fig.  3.  It  is  convenient  to  introduce  a coordinate  system 
with  origin  in  the  center  of  the  nozzle  exit  plane  and  with  the  X3  axis 
coincident  with  the  exhaust  jet  axis.  This  is  consistent  with  the  notation  of 
Section  II.  Sound  pressure  data  can  be  taken  with  either  of  two  microphone 
sweep  arms:  a rear  arm  sweeping  in  the  nozzle  exit  x^,x2  plane  and  a 
side  arm  sweeping  in  a plane  containing  the  jet  axis,  the  Xj,x3  plane.  For 
normal  operation,  the  microphones  mounted  on  the  arms  are  at  a constant  distance 
of  13  feet  from  the  center  Df  the  nozzle  exit  plane  during  a sweep.  At  this 
distance,  the  microphones  are  in  the  far  field  for  all  frequencies  greater 
than  about  100  Hz.  The  sweep  arms  are  adjustable,  so  that  the  distance  from 
the  microphone  to  the  nozzle  center  can  be  varied  over  a wide  range.  In 
particular,  the  side  arm  can  be  adjusted  to  sweep  the  microphone  in  the 
x-^Xj  plane  at  about  one  foot  from  the  coordinate  system  origin.  This  feature 
was  used  in  obtaining  near  field  data  in  the  present  study. 


3.2  ACOUSTIC  DATA 

All  acoustic  test  data  in  the  present  study  were  taken  on  high  velocity 
jets  exhausting  from  a converging  circular  nozzle  with  a 2.82  inch  exit 
diameter.  In  order  to  avoid  noise  sources  associated  with  shocks  in  the 
flow,  the  exit  Mach  number  of  the  Jet  exhaust  was  kept  below  one.  Four 
different  combinations  of  stagnation  pressures  and  temperatures  were  used 
in  order  to  vary  the  jet  exit  velocity,  and  consequently,  sound  source  power 
over  a wide  range.  The  flow  conditions  used  in  this  study  are  given  in 
the  table  below. 


TABLE  1.  TEST  CONDITIONS 


Condition  No. 

P /P 

0 1 a 

Mach  No. 

VR 

Velocity  ft/sec 

1 

1.387 

0.7 

540 

755 

2 

1.691 

0.9 

530 

943 

3 

1.691 

0.9 

1160 

1396 

4 

1.787 

0-95 

1260 

1524 
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For  each  of  the  conditions  shown  in  Table  1,  a series  of  acoustic 
data  were  obtained.  Far  field  data  in  the  x^,x3  plane  were  taken  first 
with  the  microphone  located  13  ft  from  the  origin  of  the  coordinate  system. 

One  third  octave  sound  pressure  spectra  were  obtained  at  5 degree  increments 
of  polar  angle,  <J>  , starting  at  <J)  *=  20°  to  = 90°.  The  polar  angle,  $ , 
as  shown  in  Fig.  2,  is  measured  from  the  X3  axis,  the  axis  of  the  exhaust 
jet.  Figs.  4,  5,  and  6 are  examples  of  these  far  field  spectra  measured  for 
Test  Condition  3.  From  spectrograms  such  as  these,  the  1/3  octave  far  field 
directivity  patterns  were  plotted  for  center  frequencies  of  100,  200,  500, 
1000,  2000,  5000,  and  10,000  Hz.  Figs.  7,  8,  and  9 are  examples  of  the  far 
field  directivites  for  Test  Condition  3-  It  should  be  noted  that  the  1/3 
octave  frequency  analysis  is  a constant  percentage  analysis  with  the  filter 
bandwidths  equal  to  23  percent  of  the  center  frequency.  Thus,  the  frequency 
bandwidth  of  the  100  Hz  data  in  Fig.  7 is  23  Hz  and  the  bandwidth  of  the 
1000  Hz  data  in  Fig.  8 is  230  Hz,  Data  such  as  shown  in  Figs.  7,  8,  and  9 
were  used  to  obtain  the  strengths  of  the  uncorrelated  elementary  sources, 
as  described  in  Section  II. 

Near  field  data  for  the  test  conditions  listed  in  Table  3 were  obtained 
by  readjusting  the  side  sweep  arm  to  bring  the  microphone  closer  to  the 
center  of  the  nozzle  exit  plane.  Data  were  obtained  in  the  x-^x^  plane  for 
radial  distances  between  1 ft  and  4 ft  at  approximately  1 ft  intervals.  At 
each  radial  distance,  1/3  octave  spectra  were  taken  at  several  angular 
positions  so  the  radial  variation  of  sound  pressure  could  be  plotted  for 
different  frequencies  and  angular  orientations.  Figs.  10,  11,  and  12  are 
examples  of  these  data  for  Test  Condition  3.  The  straight  solid  lines  in 
these  figures  corresponds  to  a far  field  variation  of  p^l/r2. 
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IV.  RESULTS  OF  STATIONARY  SOURCE  MODEL 


In  this  section  a comparison  of  the  results  of  the  stationary  source 
model  with  the  acoustic  data  obtained  in  the  experimental  program  are  pre- 
sented. The  first  group  of  comparisons  in  Figs.  13  through  15  illustrate  the 
degree  to  which  the  far  field  directivity  patterns  can  be  represented  by  a 
combination  of  uncorrelated  singularities  through  the  octupole  level.  The 
curves  in  these  figures  were  obtained  from  the  equation 


A^  sin 


+ A„  sin  <j>  cos  <f>  + A 


sin2  <J>  cob^*  <|) 


+ A,  cos 


i * r t 

where  A^,  Aj,  A^,  A^  are  four  positive  constants  selected  to  give  the  best  fit 
in  least  squares  sense  to  the  experimental  directivity  pattern.  The  acoustic 
pressures  in  these  and  all  other  figures  in  this  report  are  given  as  Sound 
Pressure  Levels  in  a decibel  scale. 


2 

SPL  = 10  log10  , dB 

P ref 

where 

Pref  = 2 x 10-5  N/ra2. 


Using  the  procedure  described  in  Section  II,  the  moments  of  the  largest  embedded 
monopole,  P]_,  was  determined  as 


t i « i 

P1  = min  (Alt  A2/3,  A3/3,  A^) . 


The  dipole  moments  D^,  D2 , the  quadrupole  moments  Q^,  Q,,  Q^,  and  the  reduced 
octupole  moments  A.,  A2,  A,,  A,  were  obtained  as  described  in  Section  II. 

These  values  are  given  in  the  legends  of  Figs.  13  through  15. 

These  results,  and  others  similar  to  them,  show  that  the  raultipole 
expansion  for  a stationary  source  distribution  gives  a good  representation  of  _ 
the  far  field  directivity  produced  by  a circular  Jet.  Curves  similar  to  those 
shown  in  Figs.  13  through  15  were  obtained  for  the  four  test  conditions  given 
in  Table  1 and  for  seven  frequencies:  100,  200,  500,  1000,  2000,  5000,  and 
10,000  Hz.  For  all  test  conditions  and  all  frequencies  considered,  the 
discrepancy  between  the  experimental  data  and  the  best  fit  for  the  multipole 
expansion  was  less  than  A dB. 

A comparison  of  the  radial  distribution  of  sound  pressure,  as  computed 
from  the  multipole  expansion  and  measured  in  the  experimental  program,  is 
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presented  in  Figs.  16  through  18.  Using  the  coefficients  listed  in  Figs. 

13  through  15  in  the  full  multipole  expansion,  including  the  near  field 
terms,  the  sound  pressure  level  was  computed  along  rays  emanating  from  the 
center  of  the  nozzle  exit  plane.  Examination  of  the  results  shown  in  Figs. 

16  through  18  shows  that  the  multipole  expansion  predicts  the  sound  field 
reasonably  well  at  high  frequencies  where  the  near  field  terms  are  of  minor 
importance.  As  the  frequency  under  consideration  decreases  and  the  near 
field  terms  in  the  multipole  expansion  become  increasingly  important,  the 
measured  and  predicted  sound  pressure  values  near  the  jet  exhaust  diverge. 

The  computed  values  of  the  sound  pressure  are  in  almost  all  cases  greater 
than  those  measured.  This  trend  suggests  that  the  stationary  source  model 
is  giving  too  large  a near  field  contribution  to  the  calculated  sound  field. 
Reference  to  Appendix  A shows  that  importance  of  the  near  field  terms  increases 
with  the  order  of  the  pole:  the  monopole  has  no  near  field;  the  near  field 
termB  of  the  dipole  are  proportional  to  l/(kx)2  *■  (X/2ttx)2  ; the  quadrupole 

near  field  terms  are  of  order  l/(kx)4;  and  the  octupole  terms  are  of  order 
l/(kx)6  . Thus,  the  stationary  source  model  is  resulting  in  too  great  a 
contribution  from  the  higher  order  poles.  This,  in  turn,  is  due  to  the 
directionality  of  the  far  field  sound  distribution  which  can  be  adequately 
modeled  only  with  the  inclusion  of  these  higher  order  terms. 

Closer  examination  of  the  measured  radial  sound  pressure  distribution 
for  low  frequencies  reveals  that  the  mean  square  sound  pressure  drops  off  even 
more  slowly  than  p2  03  1/r2  predicted  for  the  far  field.  This  is  perhaps  best 
seen  in  Figs.  10  and  11  where  the  straight  lines  represent  the  far  field 
relation.  This  discrepancy  obviously  cannot  be  overcome  by  the  inclusion 
of  near  field  terms,  but  must  be  accounted  for  in  some  other  way. 

In  the  next  section,  modifications  of  the  stationary  source  model  are 
introduced,  which  result  in  substantial  improvement  in  the  agreement  with  the 
experimental  data. 
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V. ' MODIFICATION  OF  STATIONARY  SOURCE  MODEL 


5.1  SOURCE  CONVECTION  EFFECTS 

In  the  stationary  source  model  presented  in  Section  II,  the  multipole 
expansion  was  used  as  if  sound  source  distribution  were  at  rest  with  respect 
to  the  observer.  In  fact,  the  elementary  sound  sources  are  fluctuating 
turbulent  eddies  which  are  being  convected  by  the  mean  flow  with  a convection 
velocity,  Vc.  This  source  motion  tends  to  beam  the  sound  field  forward  and 
produce  a directionality  that  is  not  associated  with  the  sources  themselves. 
Thus,  in  the  stationary  source  model,  the  source  synthesis  procedure  Is 
modeling  an  apparent  source  distribution  rather  than  the  true  source  distri- 
bution. 


The  jet  noise  theories  of  Lighthill  (Ref.  5)  and  Ffowcs-Williams  (Ref.  6) 
indicate  that  the  far  field  sound  intensity  radiated  from  an  axi-syiranetric 
turbulent  jet  is  given  by 


l(x,d))  = 


P V V 

-m- -J  _ 

5 2 

pax 
Ko  o 


f (j>) 


(1  - M cos  (J;)3 
c 


(5.1) 


where  x and  <p  are  the  polar  coordinates  of  the  far  field  observation  point, 
p is  a mean  jet  nixing  region  density,  V is  the  jet  exit  velocity,  and 
Mc  = Vc/aD  is  a convection  Mach  number.  Tne  convection  velocity,  Vc,  is  less 
than  the  jet  exit  velocity  and  is  thought  to  be  frequency  dependent.  The 
factor  (1  - Mc  cos  d>;  takes  into  account  the  directionality  of  the  sound 
field  caused  by  the  motion  of  the  acoustic  sources.  The  function  f(0)  is 
associated  with  the  inherent  directionality  of  the  source  distribution. 
Considering  the  failure  of  the  stationary  source  model  to  accurately  predict 
the  near  field  sound,  it  was  decided  to  use  the  multipole  expansion  procedure 
to  approximate  the  source  directionality  function  f(0). 

To  carry  out  the  revised  procedure,  the  measured  far  field  intensities 
were  divided  by  the  convection  factor  (1  - Mccos  9)-5  prior  to  determination 
of  the  multipole  moments.  This  removes  the  effect  of  source  motion  from  the 
far  field  sound,  so  that  the  multipole  expansion  includes  only  the  inherent 
directionality  of  the  source  region.  This,  in  turn,  reduces  the  weight  of 
the  higher  order  terms  in  the  expansion  with  their  strong  near  field  compon- 
ents, The  computed  sound  field  from  the  source  distribution  is  then 
multiplied  by  (1  - Mc  cos  0)  5 to  obtain  the  predicted  sound  field.  The 
unknown  in  this  procedure  is  the  ratio,  Rc,  of  the  convection  velocity  to 
the  jet  exit  velocity.  This  ratio  was  adjusted  at  each  frequency  to  remove 
as  much  angular  variation  from  the  sound  field  as  possible. 

. Figs.  19  through  22  show  some  comparisons  of  far  field  sound  directivity 
patterns  from  Test  Condition  3 with  the  results  of  the  multipole  expansion 
after  the  convection  effects  have  been  removed.  In  general,  a better  fit  to 
the  data  was  obtained  with  the  convection  factor  removed-  The  values  of 
ratio  R^  shown  in  the  legends  of  these  figures  are  not  unreasonable. 
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5.2  EXTENDED  SOURCE  EFFECTS 

In  the  stationary  source  model  presented  in  Section  II,  for  each 
frequency  the  distributed  noise  source  of  the  turbulent  jet  was  replaced 
by  an  assemblage  of  uncorrelated  elementary  sources  located  at  the  center  of 
the  nozzle  exit  plane.  However,  it  is  well  known  from  jet  noise  theory  (Ref.  7) 
that  the  location  of  the  sound  sources  is  frequency-dependent.  The  low 
frequency  sound  is  generated  by  large  scale  turbulent  eddies,  which  are 
located  further  downstream  than  the  small  scale  sources  of  high  frequency 
6ound.  In  attempting  to  account  for  this  effect*  it  was  assumed  that  the 
effective  source  for  each  1/3  octave  frequency  band  can  be  considered  as 
a point  on  the  jet  axis  and  that  the  location  of  the  source  is  primarily  a 
function  of  the  Strouhal  number,  St  = f Vj/D. 

The  Initial  investigation  of  this  effect  was  performed  by  computationally 
varying  the  source  location  and  determining  the  root  mean  square  difference  be- 
tween the  measured  sound  pressure  level  and  that  computed  from  the  multipole 
expansion  in  the  near  field.  The  coefficients  in  the  multipole  expansion  were 
determined  from  the  far  field  measurements,  after  the  convection  effects  had 
been  removed  as  described  above.  A typical  result  of  this  computation  is 
shown  in  Fig.  23.  The  approximate  source  location  xa  is  clearly  revealed 
in  this  plot  by  the  minimum  error. 

The  approximate  source  locations,  as  determined  by  the  above  procedure 
for  several  frequencies  and  flow  conditions,  are  plotted  versus  Strouhal 
Number  in  Fig.  24.  The  data  exhibited  in  this  graph  show  the  expected  trend 
that,  as  the  frequency  increases,  the  source  location  moves  closer  to  the 
nozzle  exit  plane,  x3  = 0.  The  solid  line  in  this  figure  was  drawn  by 
inspection  to  "best  fit”  the  data  and  is  given  by 


x 

s 


D 


= .7503  (St) 


-.5624 


(5.2) 


McGregor  and  Simcox  (Ref.  8 ) have  determined  peak  power  source 
location  data  for  flow  from  a circular  nozzle  by  a wall  isolation  technique. 
Their  results  show  the  same  general  trend  as  that  of  Fig.  24,  but  locate  the 
source  somewhat  farther  downstream  for  a given  Strouhal  number. 


5.3  RESULTS  OF  THE  MODIFIED  SOURCE  MODELING  PROCEDURE 

Computations  of  the  sound  field  for  Test  Conditions  1 through  4 were 
carried  out  using  the  modified  source  modeling  procedure  described  above. 

The  results  of  these  calculations  and  a comparison  with  measured  data  are 
presented  in  Figs.  25  through  44,  In  performing  these  calculations,  the 
coefficients  in  the  multipole  expansion  were  obtained  from  far  field  data  at 
x = 13  ft,  assuming  a source  location  of  x = 0.  In  all  cases  the  multipole 
expansion  was  carried  through  the  octupole  terms  and  the  convection  effects 
were  removed  before  evaluation  of  the  multipole  moments.  In  Figs.  25 
through  44  all  angles  and  distances  are  referred  to  a coordinate  system 
fixed  in  the  nozzle  exit  plane,  as  shown  In  Fig. 2,  although  the  calculations 
using  the  equations  of  Section  II  were  made  with  angles  and  distances  based 
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on  the  source  location  for  each  frequency.  No  comparisons  are  presented 
for  frequencies  above  2000  Hz,  since  at  higher  frequencies  the  sources  are 
located  very  near  the  no2zle  exit  and  all  microphone  positions  are  in  the 
far  field. 

In  general,  the  predicted  values  agree  reasonably  well  with  the 
microphone  data.  However,  fairly  large  discrepancies  still  occur  for 
measurements  very  close  to  the  nozzle  and  for  small  angles  to  the  jet  axis. 
The  microphone  measurements  at  90°  are  those  which  are  best  described  by 
the  prediction  scheme.  This  is,  perhaps,  to  be  expected,  Bince  in  this 
direction  the  radiated  sound  is  least  affected  by  source  convection  or 
flow  interaction  effects. 
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VI.  SUMMARY  AND  CONCLUSIONS 


At  the  outset  of  this  investigation,  it  had  been  hoped  that  the  acoustic 
field  produced  by  a free  circular  jet  could  be  well  modeled  by  a stationary 
distributed  source,  which  could  in  turn  be  closely  approximated  by  a multipole 
expansion  centered  at  the  origin  of  the  coordinate  system.  It  turned  out  that 
this  approach  produced  a near  field  contribution  much  larger  than  obtained  by 
actual  measurements.  In  order  to  achieve  better  agreement  with  the  experimen- 
tal results,  it  was  necessary  to  use  a more  realistic  model  for  the  source 
region.  In  particular,  the  fact  that  the  acoustic  sources  in  a jet  flow  are 
moving,  being  convected  by  the  mean  flow,  had  to  be  taken  into  account.  Also, 
the  fact  that  the  low  frequency  noise  sources  are  located  further  downstream 
from  the  nozzle  exit  plane  than  the  high  frequency  sources  had  to  be  incor- 
porated into  the  jet  noise  source  model.  With  these  modifications,  calculations 
(of  the  near  acoustic  field  were  carried  out.  A direct  comparison  of  the 
calculated  and  measured  sound  field  was  presented  in  Section  V and  justify  the 
following  conclusion:  The  near  and  far  acoustic  fields  produced  by  a circular 
jet  can  be  reasonably  well  calculated,  using  a model  source  region  consisting 
of  uncorrelated  elementary  sources  if  the  effects  of  source  motion  and 
distribution  with  frequency  are  taken  Into  account. 

A second  conclusion  reached  on  the  basis  of  this  investigation  is  that 
predicting  free  field  acoustic  data  from  measurements  in  the  test  cell  will  be 
very  difficult.  This  conclusion  is  in  part  based  on  the  fact  that  there 
appears  to  be  no  very  strong  direct  near  field  contribution  that  is  simply 
related  to  the  far  acoustic  field.  This  being  the  case,  it  will  be  very 
difficult  to  separate  the  direct  radiation  field  from  the  reverberant  field 
caused  by  reflections  from  the  test  section  wall. 

A further  complication  in  the  test  cell  is  the  fact  that  the  engine 
exhaust  is  captured  by  a diffuser  or  scavenger  duct  a short  distance  downstream 
of  the  nozzle  exit.  Thus,  much  of  the  sound  generated  by  the  exhaust  jet  may 
not  be  escaping  into  the  test  cell.  On  the  other  hand,  if  the  nozzle  exhaust 
impinges  or  scrubs  the  diffuser,  additional  noise  sources  will  be  present. 

There  can  be  no  doubt  that  the  sound  field  produced  by  the  exhaust  jet  - diffuser 
combination  is  quite  different  from  that  produced  by  the  jet  alone. 

It  is  thought  that  a better  understanding  of  the  acoustic  characteristics 
of  the  test  cell  and  an  investigation  of  the  effect  of  the  diffuser  on  the  jet 

sound  field  must  precede  further  efforts  to  predict  far  field  acoustic  data  from 

measurements  in  the  test  cell.  A study  of  these  factors  would  also  be  of 
considerable  value  in  any  effort  to  predict  or  assess  the  acoustic  environment 
in  a test  cell.  For  these  reasons  it  Is  recommended  that  Phase  II  of  the  present 
contract  be  directed  toward  a detailed  investigation  of  the  acoustic  field 

produced  by  the  interaction  of  the  jet  exhaust  with  the  scavenger  duct  and  an 

attempt  to  theoretically  describe  the  sound  field  within  a test  cell. 
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APPENDIX  A 


The  mean  square  acoustic  pressures  associated  with  the  monopole,  dipole, 
quadruopole  and  octupole  type  sources  are  given  below: 


Monopole: 


-5-  M2 


000 


(4TTX)' 


(no  nearfield  terms) 


Dipole: 


p2.  - — ^—5-  [1  + — ^2  1 <m2ioo  s1r2  + + m2ooi  cos2  W 

d (4*x)2  (kx)2  100  001 


Quadrupole: 


kv 

2 n 200  4 1 4 2 

p - =*5-  [sin4  if  + — ^ (3  sin  if  - 4 sin  $ + 2) 

q (4TT*r  (kx)z 


+ 1 , (9  sin4  $ - 6 sin2  <f  + 2)  ] 

(kx)4 


kV, 


+ =-=-  [cos  if  + r-  (3  cos  $ - 4 cos  if  + 1) 
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— — — (9  cos  if  - 6 cos  if  + 1) 
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kV 
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Octupole: 


k6M2 

p20  «=  [sin6  4 + — ^ (sin6  4-3  sin4  4+2  sin2*) 


(4irx)‘ 


(kx)‘ 


+ — 2 ■ (sin6  4-72  sin4  4+36  sin2  4) 

(kx)4 
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APPENDIX  B 


Figures 
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Fig.  4.  Far  Field  Spectrum  for  Test  Condition  3 at  <j>  “ 30". 
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Fig.  6.  Far  Field  Spectrum  for  Test  Condition  3 at  t = 90°. 
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One  Third  Octave  Far  Field  Directivity  for  a Center  Frequency  of  100  Hz.  Test  Condition  3 
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Fig.  8.  One  Third  Octave  Far  Field  Directivity  for  a Center  Frequency  of  1000  Hz.  Test  Condition  3 
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One  Third  Octave  Far  Field  Directivity  for  a Center  Frequency  of  10,000  Hz.  Test  Condition  3. 
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Fig.  10.  Radial  Variation  of  Sound  Pressure  Level  for  Test  Condition  3 at  100  Hz. 
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Fig,  11.  Radial  Variation  of  Sound  Pressure  Level  for  Test  Condition  3 at 
1000  Hz ■ 
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Fig.  13.  Comparison  of  Measured  Sound  Pressure  Level  Directivity  with  Least  Squares  Fit  for  Test 
Condition.  2 at  100  Hz. 
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Fig.  14.  Comparison  of  Measured  Sound  Pressure  Level  Directivity  with  Least  Squares  Fit  for  lest 
Condition  2 at  1000  Hz. 
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Fig.  15.  Comparison  of  Measured  Sound  Pressure  Level  Directivity  with  Least  Squares  Fit  for  Test 
Condition  2 at  10,000  Hz. 
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Fig.  16.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  for 
Test  Condition  2 at  100  Hz. 
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Fig.  17.  Comparison  of  Measured  and  Predicted 
Test  Condition  2 at  1000  Hz. 
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Fig.  18.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels 
for  Te.st  Condition  2 at  10,000  Hz. 
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Fig.  19.  Comparison  of  Measured  Far  Field  Sound  Directivity  after  Removal  of  Convection  Effects  with 
Least  Squares  Fit  for  Test  Condition  3 at  100  Hz. 
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Fig.  20.  Comparison  of  Measured  Far  Field  Sound  Directivity  after  Removal  of  Convection  Effects  with 
Least  Squares  Fit  for  Test  Condition  3 at  200  Hz. 
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Fig.  21.  Comparison  of  Measured  Far  Field  Sound  Directivity  after  Removal  of  Convection  Effects  with  Least 
Squares  Fit  for  Test  Condition  3 at  500  Hz. 
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Fig.  22.  Comparison  of  Measured  Far  Field  Sound  Directivity  after  Removal  of  Convection  Effects  with 
Least  Squares  Fit  for  Test  Condition  3 at  1000  Hz. 


AEDC-TR-79-3B 


Test 

Condition 

Mach  No. 

Stagnation  Tamp  (*F) 

□ 

4 

0.95 

BOO.O 

A 

1 

0.70 

80.0 

O 

3 

0.90 

700.0 

Fig.  24.  Optimum  Source  Location  aa  a Function  of  Strouhal  Number 


AEDC-TR  79  36 

FREQUENCY [HZ!*  100 

Q 

40.  DEG. 

T = 80  DEG  F 

0 

0 

60.  DEG. 

HRCH  NO.*  .70 

A 

90.  OEG. 

— 

PREDICTED  SPL  AT 

INDICATED  ANGLES 

2.99  3.00  4.00  6.00  B.00  20.00  25.00 

DISTANCE  FROM  NOZZLE, FT 


Fig.  25.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  1 at  100  Hz 
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Fig.  29.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 

Modified  Source  Modeling  Technique  for  Test  Condition  1 at  2000  Hz. 
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Fig.  30.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  2 at  100  Hz. 
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Fig.  31.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  2 at  200  Hz. 
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Fig.  34.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 

Modified  Source  Modeling  Technique  for  Test  Condition  2 at  2000  Hz. 
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Fig.  37.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  3 at  500  Hz. 
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Fig.  38.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 

Modified  Source  Modeling  Technique  for  Test  Condition  3 at  1000  Hz. 
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Fig.  39.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 

Modified  Source  Modeling  Technique  for  Test  Condition  3 at  2000  Hz 
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Pig.  40.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  4 at  100  Hz. 
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Fig.  41.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  4 at  200  Hz. 
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Tig.  42.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  4 at  500  Hz. 
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Fig.  43.  Comparison  of  Measured  and  Predicted  Sound  Pressure  Levels  using 
Modified  Source  Modeling  Technique  for  Test  Condition  4 at  100  Hz. 
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LIST  OF  SYMBOLS 

6 2 2 

Octupole  coefficients,  k /M  » etc. 

Ambient  speed  of  sound 

Dipole  coefficient,  kVio0/(W  ’,  etc. 

Nozzle  diameter 
Source  function 
Frequency 

Source  directivity  function 
Imaginary  unit,  i/^T 

~ 

Far  Field  intensity,  p /pa^ 

Wave  number,  to/a^ 

Convection  Mach  number,  V^/ 

Multipole  moment 

Multipole  moment  amplitude 

2 2 

Monopole  coefficient,  M qqq/C%x) 

Stagnation  pressure 
Ambient  pressure 
Acoustic  pressure 

Acoustic  pressure  associated  with  monopole  terms 

Acoustic  pressure  associated  with  dipole  terms 

Acoustic  pressure  associated  with  quadrupole  terms 

Acoustic  pressure  associated  with  octupole  terms 

Quadrupole  coefficients,  k4M22oo/(4ll302  etc. 

Distance  from  source  point  to  observation  point,  |x  - y| 

Convection  velocity  ratio, 

Radius  In  cylindrical  coordinates 

~~2  2 

Sound  pressure  level,  10  log  p /p 
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AEDC-TR-79-36 


St 


x1,x2,x3 


yl,y2,y3 

6 

X 

p 

u 


Strouhal  number,  fV^/D 
Stagnation  temperature 
Convection  velocity 
Jet  exit  velocity 
Coordinates  of  field  point 
Source  location 
Distance  to  field  point 
Coordinates  of  source  point 
Azimuthal  angle 
Wavelength,  a^/f 
Dens ity 
Polar  angle 
Frequency,  2irf 
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